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Based on a thorough numerical analysis of the spectrum 
of Elarper's operator, which describes, e.g., an electron on a 
two-dimensional lattice subjected to a magnetic field perpen- 
dicular to the lattice plane, we make the following conjecture: 
For any value of the incommensurability parameter a of the 
operator its spectrum can be covered by the bands of the 
spectrum for every rational approximant of g after stretching 
them by factors with a common upper bound. We show that 
this conjecture has the following important consequences: For 
all irrational values of a the spectrum is (i) a zero measure 
Cantor set and has (ii) a Ifausdorff dimension less or equal 
to 1/2. We propose that our numerical approach may be a 
guide in finding a rigorous proof of these results. 

PACS numbers: 73.20.Dx, 71.30.+h, 02.30. Tb 



I. INTRODUCTION 

The Harper model [Q may be considered as the most 
important system in the realm of quasi-periodic sys- 
tems. It is described by the one-dimensional tight- 
binding Hamiltonian 

Hx^a.u XI "«+i°^n + + Acos(27rcrn + v)aia^, 

n 

(1) 

where A, cr, and v are real valued parameters, and ajj and 
a„ are creation and annihilation operators at site n. De- 
pending on whether a is rational or irrational this Hamil- 
tonian is periodic or quasi-periodic, resp. The model ex- 
hibits rich spectral features which showed up for the first 
time in the numerical work of Hofstadter Q| who com- 
puted its spectrum for A = 2 as a function of a. The 
beautiful resulting graph (Fig. 1) is today known as Hof- 
stadter's butterfly and has attracted a lot of attention 
among physicists and mathematicians ever since. 

From a physical point of view the Harper model was 
first introduced to describe electrons moving in a two- 
dimensional periodic potential with a square unit cell and 
subjected to a weak magnetic field perpendicular to the 
potential plane In this case. Harper's model de- 

scribes the eigenfunctions in one direction of the under- 
lying Bravais lattice after separating plane waves in the 
perpendicular direction. Then, the parameter A is twice 
the ratio of the modulation amplitudes in these two direc- 
tions, a is given by the number of magnetic flux quanta 



per unit cell of the potential, and v is the wave num- 
ber of the plane wave being separated. Such a system 
can be realized, e.g., by lateral superlattices on semicon- 
ductor heterostructures in which the first experimental 
indications of Hofstadter's butterfly were found recently 
Harper's model also appeared in the opposite limit, 
namely for electrons in a two-dimensional periodic po- 
tential with strong magnetic field Q , as well as in other 
physical contexts, e.g., phonons in a quasi-periodic po- 
tential and superconducting networks jsj . Apart from 
describing real physical systems it is an important and 
simple model exhibiting a metal- insulator transition 
if (J is a typical irrational, then for A < 2 all eigenfunc- 
tions are extended while for A > 2 they are localized. 
Hence, at A = 2 a metal-insulator transition occurs and 
unusual properties of the spectrum and the eigenfunc- 
tions are expected. This has led mathematicians to study 
Harper's model as a simple example of a system with a 
possibly singular continuous spectrum. In particular, at 
A = 2 the spectrum has been conjectured to be a zero 
measure Cantor set for irrational values of a . This 
has been shown for a large class of irrationals by Helf- 
fer and Sjostrand ||ll| and Last ||l^, which, however, do 
not cover the irrationals typically used in numerical in- 
vestigations, namely irrationals with bounded continued 
fraction expansions, e.g., the Golden Mean (V5-1) /2. So 
still the question remains whether the spectrum is a zero 
measure Cantor set for all irrational values of a . 

For both, physicists and mathematicians, the Harper 
model is also an important example for the study of the 
relation between the spectral properties and the dynam- 
ics of quantum systems |]l5|-^l|. In 1988, Hiramoto and 
Abe numerically discovered anomalous diffusive spread- 
ing of wave packets in the Harper model for A = 2 [p^ . 
Furthermore, the model has been used for testing rela- 
tions proposed between the fractal properties of the spec- 
tral measure and the temporal decay of autocorrelations 
]l6t as well as for the asymptotic growth of the moments 
of a wave packet ||l^,|2^,^. In these studies one is thus 
interested in fractal properties of the spectrum, the eas- 
iest being the Hausdorff dimension. First, it was investi- 
gated numerically | p^ , p^ -p^ and recently a perturbative 
calculation has been presented for a the Golden and the 
Silver Mean |2^ . For a class of irrationals, which has zero 
Lebesgue measure, Last proved the Hausdorff dimension 
to be bounded from above by 1/2 |jl^. It remains open 
whether this bound is valid for all irrationals. 
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In this article we study the operator at the critical 
point A = 2. Based on a numerical analysis we present 
our conjecture that there is a constant R such that the 
spectrum for any value of a can be covered by the q bands 
of the spectrum for any rational approximant p/qoiu al- 
ter stretching each of these bands by a factor smaller than 
R. This conjecture permits us to prove for all irrational 
a that the spectrum is (i) a zero measure Cantor set and 
has (ii) a Hausdorff dimension less or equal to 1/2. After 
presenting the relevant basic facts on Harper's operator 
in Sec. II, in Sec. Ill we precisely state our conjecture and 
deduce the above mentioned statements on the measure 
of the spectrum and its Hausdorff dimension. Section IV 
contains a thorough numerical analysis of Hofstadter's 
butterfly substantiating the validity of our conjecture. 

II. BASIC FACTS 

In the following we review the facts about the spectrum 
of Harper's operator [Eq. (^] that are used in the proofs 
presented in Sec. III. The relevant energy spectrum S{a) 
of the Harper model for A = 2 and a fixed value of a is 
given by the set 



S{<y) Spec {H2.a,y) , 



(2) 



where Spec{H2^a,u) denotes the spectrum of i?2,cr.i/- Fig- 
ure 1 shows S{(t) vs. ct, the so-called Hofstadter butter- 
fly revealing its symmetries with respect to i? = 
and a = 1/2. 

For rational a = p/ q, with p and q relatively prime, 
the operator (|^) is periodic and 5(<t) consists of q bands. 
These bands are separated hy q~ \ gaps in the case of q 
odd and q—2 gaps for q even, for which the middle bands 
touch each other a.t E = The Lebesgue measure 
\S{p/q)\ of the spectrum obeys [|l2[ 



2(75-1-1) 



< 



< 



(3) 



where e — exp(l) = 2, 71 . . . 

For irrational values of a, the operator is quasi- 
periodic and its spectrum contains no isolated points 
as follows from general results on ergodic Jacobi matri- 
ces [^. Furthermore, it is independent of ||2^, such 
that S{cr) is a closed set. Thus, if S{a) has measure zero 
it will be nowhere dense and therefore will be a Cantor 
set. As the measure \S{p/q)\ of rational approximants 
p/q of a decreases to zero with increasing q [Eq. (H)], one 
may naively think that indeed |>5'(f7)| — holds. Rig- 
orously, however, |S'((j)| can only be shown to vanish for 
two classes of irrationals a: One class consists of the irra- 
tionals for which all coefficients in the continued fraction 
expansion are larger than some constant pT[ | , while the ir- 
rationals in the other class have an unbounded continued 



fraction expansion |l|. The latter proof uses a continu- 
ity property of the spectrum, namely that the spectrum 
S{a) is Holder continuous of order 1/2 i.e., for ev- 
ery eigenvalue E in the spectrum for a there is a nearby 
eigenvalue E' in the spectrum for a' with 

\E-E'\<6{2\cr^a'\y/^. (4) 

This allows to construct a sequence of covers of S{a) us- 
ing the approximants {Pn/Qn}^=o '^^ ^^"^ continued frac- 
tion expansion of a. For each n one attaches intervals of 
length 6(2|(T — Pn/qn\y^'^ to each of the 2qn band edges 
of S{pn/qn)- According to Eq. (|^) this defines a cover 
Sq^^ (a) of S{a) consisting of qn intervals. The measure of 
any of these covers is 



I2q„ 
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1/2 



(5) 



(6) 



by Eq. (0). For n ^ oo the measure vanishes only if 
lim.n^aoq^W ~ Pn/qn\ — 0, which is fulfilled for irra- 
tionals with an unbounded continued fraction expansion 
and proves a zero measure Cantor set spectrum for these 
irrationals |12|. For general a, which fulfill 

C 



< 



q- 



(7) 



for some constant C > 0, covers as described above do 
not permit to show |S'(tT)| = 0. With respect to the Haus- 
dorff dimension these covers have been used by Last |l^] 
to show an upper bound of 1 /2 for the zero measure set of 
irrationals fulfilling qn\cF — Pn/ qn\ < C for all n and some 
constant C > 0. In the following, we will present a con- 
jecture, stating that there is a useful cover obtained by 
stretching the bands of a rational approximant of a. As- 
suming this conjecture we prove that the Lebesgue mea- 
sure of the spectrum is zero as well as that its Hausdorff 
dimension is less or equal to 1/2 for all irrationals. 



III. THE CONJECTURE AND ITS 
IMPLICATIONS 

In this section we will show how the zero measure Can- 
tor structure of spectra for irrational values of a and an 
upper bound of their Hausdorff dimension follow from 
our 

Conjecture 1 : There is a constant R > such that 
for any irrational cr G [0,1] and any two rational approx- 
imants p/q and p' /q' with q < q' , which are obtained by 
truncating the continued fraction expansion of a , there 
exists a cover Sq{p' / q') of S{p' , q') consisting of q inter- 
vals with the property 



< R 



S 



(8) 
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Before giving the explicit construction of such a cover and 
numerically substantiating this Conjecture in the next 
section, we will now discuss its consequences. One can 
replace p'/g' in Eq. (^) by an irrational a as is stated in 
Corollary 1 : Assuming Conjecture 1, there is a con- 
stant R > such that for any irrational a G [0, 1] and any 
rational approximant p / q obtained by truncating the con- 
tinued fraction expansion of a there exists a cover Sq (a) 
of S{a) consisting of q intervals with the property 



< 



ieR+1 

q 



(9) 



Proof. Let a be an irrational and let p/q be a ratio- 
nal approximant of a, which is obtained by truncating 
its continued fraction expansion. Then, one can always 
choose a sufficcntly high rational approximant p' /q' of 
a, which again is obtained by truncating its continued 
fraction expansion, such that 



q' > 12q'^V2C, 



(10) 



where C is the constant of Eq. (^. We now define a 
q-cover Sq{a) of the spectrum for a using the continuity 
property Eq. (^). To this end we take the q intervals of 
the cover Sq(p' /q'), as given by Conjecture 1, and enlarge 
them by 6(2|cr— (p'/q')l)^^^ the lower and upper ends, 
respectively. For the Lebesgue measure |S'g(CT)| of Sq{a) 
we thus have 



< 



< 



< 



2g • 6 2 
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(11) 

(12) 
(13) 



where we have used Eq. (|^). Finally, using Conjecture 1 
and Eq. (|) we obtain Eq. (|). □ 

From this Corollary follows 
Theorem 1 : Assuming Conjecture 1, for any irrational 
a € [0, 1] the set S{a) is of zero Lebesgue measure. Proof. 
By taking the limit g — > cxd in Eq. (^ this follows imme- 
diately. □ 

The Cantor structure of the spectrum for irrational a 
is a direct consequence of Theorem 1 and stated in 
Corollary 2 : Assuming Conjecture 1, for any irrational 
a e [0, 1] the set S{a) is a Cantor set. Proof. Since 
the spectrum S{a) is closed, Theorem 1 shows that it is 
nowhere dense. As it is also known to have no isolated 
points it follows that S{a) is a Cantor set. □ 

As a second application of Conjecture 1, we now give 
an upper bound on the Hausdorff dimension Dh of the 
spectrum. Most of the proof is contained in 
Lemma 2 (Ref. [0): Let S' C M, and suppose that S 
has a sequence of covers: {<S'g„}^i, S C Sq^, such that 



each Sq^ is a union of qn intervals, (7„ — > cx) as n —s- oo, 
and for each n: 



qn' 

where /? and c are positive constants; then: 

1 



DniS) < 



(14) 



(15) 



Proof See Ref. p 



Theorem 2 : Assuming Conjecture 1, for any irra- 
tional CT e [0,1] we have 



Dh{S{<t)) < 



1 



(16) 



Proof. Choosing for Sq^ in Lemma 1 a cover satisfying 
Eq.(^ of Corollary 1 and then comparing Eqs. (^) and 
( p^ ) yields j3 = 1. Inserting this into Eq. ( p^ ) concludes 
the proof of Theorem 2. □ 



IV. THE COVER 

In this section we present a specific cover which is nu- 
merically shown to satisfy Conjecture 1. A good start- 
ing point for covering the spectrum of Harper's operator 
[Eq. (|l|)] for some rational a' are the bands of the spec- 
trum for a rational a obtained by truncating the contin- 
ued fraction expansion of a' . In fact, there is a natural 
assignment of any of these bands to a cluster of bands of 
the spectrum for cr', which is determined by Hofstadter's 
Rules (see Appendix). Very often a band is too narrow to 
cover its related cluster, e.g., the lowest band of a = 1/3 
does not cover the related cluster consisting of the lowest 
two bands for a' — 2/5. Therefore, we introduce the fac- 
tor i?i(|, y) by which the z-th band, i = 1, 2, . . . , of 
the spectrum for a = p/ q has to be stretched in order to 
cover its related cluster of the spectrum for a' = p' / q' ■ 

We numerically analyzed the stretching factors Ri for 
all allowed pairs p/q and p' /q' with q' < 300 and found 
that they are bounded from above by 



max 

g'<300 



Ri 



p p_ 

1 I 

q q 



< 2.1 



(17) 



If a bound would exist without the restriction q' < 300 
for any two rational approximants p/q and p'/q' with 
q < q' , that are obtained by truncating the continued 
fraction expansion of some irrational tr S [0, 1], then the 
cover described above would fulfill Conjecture 1. In the 
remaining part of this section we will argue with the help 
of further numerical analysis that this is very likely to be 
the case. 
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We first investigate for which bands and field values 
the maximum stretching factor for a fixed q occurs. We 
observe that the maximum stretching factor always oc- 
curs when covering the spectrum for 1/(9 + 1) with the 
bands of the spectrum for At these magnetic fields 
the maximum stretching factor stems from the central 
bands, namely 



max 

q'<300 



p p_ 
'9' 



= R 



1 



(18) 



where [x\ denotes the integer part of x. For example, in 
the case q — S {q — A) the maximum stretching factor 
occurs when using the second band oi p/q — 1/3 {p/q = 
1 /4) to cover the cluster containing the second and third 
bandofp'/g' = 1/4 {p' /q' = 1/5) (Fig. 1). Figure 2 shows 

^l^a+ij ^ function of q. It has a maximum 

value of about 2.015 at q ~ 300 and is asymptotically 
decreasing towards 2, so that it is bounded from above. 
This behaviour is substantiated (Fig. 2) by using results 
from Thouless and Tan |29[| giving an approximation of 
the band edges for magnetic fields of the form 1/q. 

Since Hofstadter's butterfly is self-similar, one expects 
that the stretching factors coming from smaller distorted 
copies of the butterfly play an important role. In fact, 
one may wonder if they require larger stretching factors, 
because of their distortion. For example, such a copy lives 
between the lowest two bands of 1/3 and 1/2 (Fig. 1) with 
an effective magnetic field CTeff given by cr = 1/(2 4- cTcff)- 
Assuming that in the smaller copies the same bands as 
in the original butterfly need the largest stretching fac- 
tors, we investigated the stretching factors of the center 
bands of this distorted copy for cTeff = 1/ qeS when cover- 
ing the spectrum for a'^f^ = l/((7eff + !)■ As can be seen 
in Fig. 2 (crosses) they converge towards the correspond- 
ing stretching factors of the original butterfly, however, 
with deviations for small gcff- A systematic analysis of 
these deviations also for other small copies of the but- 
terfly shows that these deviations are bounded (Fig. 3). 
This suggests that Eq. ( p7[ ) remains valid even without 
the restriction q' < 300 and thereby gives strong numer- 
ical evidence for the validity of Conjecture 1. 



V. CONCLUSION 

We have conjectured that the spectrum of the Harper 
operator for a rational value of the incommensurability 
parameter a can be covered by the spectrum for a ra- 
tional approximant of a after stretching the bands of 
the latter by a factor R independent of a. This factor 
is numerically found to be R w 2.1. We showed that 
from this Conjecture follows that the spectrum is a zero 
measure Cantor set for all irrational a. Furthermore, it 
implies that the Hausdorff dimension of these spectra is 



bounded from above by 1/2. While these results are for 
irrational values of a, the underlying conjecture is stated 
for pairs of rational values of a. This allowed a detailed 
numerical study giving strong numerical evidence for our 
Conjecture. As we give an explicit construction of the 
cover used for the conjecture, we hope that our numeri- 
cal approach is helpful in finding a rigorous proof of these 
results. 
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APPENDIX 

Here we describe in detail how the cover used in Section 
IV is constructed using Hofstadter's rules (|,^. This set 
of rules describes qualitatively the self-similar structure 
of the graph in Fig. 1. For a given magnetic field a the 
energy spectrum of the Harper operator consists of three 
non-overlapping parts, one central cluster C and two side 
clusters L and R. These three parts are distorted Harper 
spectra for effective field values ctl, ac and an which are 
given by 



crc = 



■{i} forO<a<i 
.{t^} fori<.<l 



(19) 



(20) 



for ^ < cr < 1 



where {x} denotes the fractional part of x. These rules 
also apply to the distorted spectra such that the spec- 
trum splits into three clusters, each of which splits into 
three clusters. For irrational a this splitting is contin- 
ued ad infinitum. In contrast, for rational a this process 
eventually stops for each cluster when it consists of either 
a single band (cTcff = 0) or of two bands (dcff — 1/2); 
see Fig. 4 for a = 5/8. For a', where a was obtained 
by truncating the continued fraction expansion of a' , we 
simultaneously apply Hofstadter's rules up to the same 
level as for a; see Fig. 4 for cr' — 34/55. This defines a 
straightforward assignment of a band with a^s = in the 
spectrum for cr to a cluster of bands in the spectrum for 
a'. In the case cTcs — 1/2 one has two bands that have to 
cover the corresponding cluster, which consists of three 
subclusters. We assign one band to two subclusters and 
the other band to the third subcluster (Fig. 4) such that 
the possibly necessary stretching factors for the bands 
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in order to cover their subcluster are minimized. By this 
procedure we assign every band of cr to a cluster of bands 
of cr', where the union of these clusters is the complete 
spectrum for a' . 
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FIG. 1. The energy spectrum of Harper's operator [Eq.(jl|)] 
for A = 2 as a function of the incommensurability parameter 
(T. The resulting graph is known as Hofstadter's butterfly. 
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FIG. 2. Stretching factor gTi} (dia- 

monds) compared to its analytic value (solid line) obtained 
from an approximation of the band edges for magnetic fields 
of the form 1/q by Thouless and Tan [^. The crosses show 
the corresponding stretching factors of the distorted copy of 
the butterfly between the lowest two bands of 1/3 and 1/2 
(Fig. 1) vs. gefi (see text). 

FIG. 3. Maximum deviation A(g) of stretching factors in 
the copies of the butterfly around p/q from the correspond- 
ing stretching factors in the original butterfly (diamonds in 
Fig. 2). Only stretching factors for l/^cff, analogous to 
(^qi 5Tt) original butterfly, are taken into ac- 

count. For different geff the maximum is taken over all copies 
of the butterfly and all p. For each geff the deviations remain 
bounded and the upper bounds decrease with increasing q^s. 

FIG. 4. Assignment of bands in the spectrum for cr = 5/8 
to clusters of bands in the spectrum for a' — 34/55 by si- 
multaneously applying Hofstadter's rules. The assignment 
is straightforward, except for the case a^B = 1/2 where two 
bands have to be assigned to three subclusters (see text) . The 
arrows mark the cases where bands have to be stretched in 
order to cover their assigned cluster. 
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